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In this paper, we demonstrate, first in literature known to us, that potential functions can be constructed in 
continuous dissipative chaotic systems and can be used to reveal their dynamical properties. To attain this 
aim, a Lorenz-likc system is proposed and rigorously proved chaotic for exemplified analysis. We explicitly 
construct a potential function monotonically decreasing along the system's dynamics, revealing the structure 
of the chaotic strange attractor. The potential function can have different forms of construction. We also 
decompose the dynamical system to explain for the different origins of chaotic attractor and strange attrac- 
tor. Consequently, reasons for the existence of both chaotic nonstrange attractors and nonchaotic strange 
attractors are clearly discussed within current decomposition framework. 



Potential function (also known as energy func- 
tion, generalized Hamiltonian, or Lyapunov func- 
tion, under different contexts) describes nonlin- 
ear dynamical system from a global point of view. 
Along this scalar function, all the states in phase 
space move downward. Consequently, the poten- 
tial function accounts for both detailed structure 
and long term trend of the system's dynamics, in- 
dicating its performance and stability at the same 
time. It is also anticipated that when the poten- 
tial function of a system becomes constant, the 
system has evolved into an "attractor". Hence, 
potential function has special theoretical impor- 
tance to chaotic system in that it helps reveal the 
complex structure of chaotic attractor. Unfortu- 
nately, fundamental difficulties pertaining to its 
construction in nonlinear dynamical systems still 
need solution. Failure in its construction has even 
prompted claims that potential function does not 
exist in complex systems 1,2 . In this paper, we 
demonstrate that potential functions can be con- 
structed in chaotic dynamical systems. We first 
present a simplified geometrical Lorenz attrac- 
tor, and rigorously prove that it is chaotic by its 
Poincare map (which is itself an effort interest- 
ing to many researchers 3 -). Then we analytically 
construct a potential function for the system, ac- 
counting for the structure of the chaotic attrac- 
tor. With this potential function, we discover 
that chaotic attractor may not be a strange at- 
tractor and vice versa. This corresponds to pre- 
vious observations and is explained in detail by 
virtue of our constructive approach. 
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I. INTRODUCTION 

Nonlinear dynamics underlying many natural and 
technological systems is described by a set of ordinary 
differential equations: 



This description defines a fixed rule dictating the trend 
of evolution that a current state will follow into its im- 
mediate future. Overtime, systems with even simple de- 
terministic rules can give rise to seemingly "fortuitous" 
phenomenal, generally known as chaos. It remains an 
intriguing problem as to how these phenomena can be 
analyzed in a global sensed, under the light of a generic 
framework in nonequilibrium dynamics^. Ideally, such a 
generic framework should account for: a unified (proba- 
bly geometric) structure underlying the equations of evo- 
lution, a comparable measure of the system's different 
states, and an accurate reflection of the dynamic process 
generated by the system. 

On an historical account, active search for such a 
generic framework begins in the 1970s, when Rene Thorn 
and Christopher Zeeman proposed that potential func- 
tions and their deformation^ can describe "the evolution 
of form in all aspects of nature, and hence it embodies 
a theory of great generality"—. Thorn's proposal is inge- 
nious. Because such a scalar function incorporates the 
previously mentioned attributes into one single quantity: 
potential function. It not only generalizes existing ap- 
proaches of Lyapunov function and first integral, but also 
encompasses concepts like stability^ and reversibility^ 
into a uniform framework. Even, natural and technologi- 
cal systems can directly be modeled with potential func- 
tions so that the behaviors can still be described "even 
when all the internal parameters describing the system 
are not explicitly known"—. 

However, according to Stephen Smale, Thorn's math- 
ematical approach "deals in with only a few known 
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examples"—, and hence lacks practical effectiveness. 
What Thorn and his followers failed to obtain is a proof 
of existence or even a method of construction for poten- 
tial functions in complex systems, such as systems with 
oscillation or chaotic behaviors. Frustration has even 
prompted claims that potential functions do not exist 
in the situations of complex dynamical systems^ 

Previously, we have rigorously defined potential func- 
tion in mathematical terms (see definition [T]) and al- 
ready demonstrated that potential functions (or Lya- 
punov functions) can be analytically constructed in os- 
cillating system a 12 ! 13 . In this paper, we further motivate 
such research by showing the construction of potential 
functions in chaotic systems, and providing additional 
insights for chaotic and strange attractors. 

Actually, constructing potential-like functions in 
chaotic systems, functions with a restricted part of the 
properties held by potential functions, is an approach al- 
ready taken by researchers. Until very recently, there 
are still various efforts addressing the issue. There are 
generalized Hamiltonian approach^, energy-like func- 
tion technique^, minimum action method^, and etc., 
in search for a unified description of chaotic dynam- 
ics. These previous methods all construct a potential- 
like function to analyze some chaotic system such as the 
Lorenz system^. Unfortunately, the scalar functions in 
these works all lack certain important properties (see sec- 
tion 6). 

Those important drawbacks of the existing methods 
motivate a real potential function to describe the behav- 
ior of chaotic systems. Ideally, it should be a continu- 
ous function in the phase space, monotonically decreasing 
with time. Also, when time approaches infinity, poten- 
tial function should stabilize to some finite quantity if 
the original system is not divergent. In this paper, we 
combine these ideas together as a potential function and 
try to obtain it in chaotic systems. 

The paper is organized as follows. First of all, we 
define an ideal potential function for the description of 
chaotic systems. To analyze chaotic systems in detail 
with this potential function, we create an attractor that 
is chaotic by definition. Then, potential function for this 
chaotic attractor is constructed, showing the structure of 
the chaotic strange attractor. In addition, our framework 
provides a decomposition of the original vector field. The 
decomposition helps understand the different origins for 
chaotic attractor and strange attractor, explaining why 
there exists both chaotic nonstrange attractors and non- 
chaotic strange attractors. 



II. POTENTIAL FUNCTION 

We first state the definition of a potential function, 
which is a natural description of dynamical systems with 
monotonic properties. Then we will discuss a decomposi- 
tion scheme of generic dynamical systems associated with 
the potential function. 



Definition 1 (Potential Functioni&±2) . Let * : R" -> R 

be a continuous junction. Then ^ satisfying the following 
condition is called a potential function for the dynamical 
system x = f(x) : R" ->■ R™. 

(a) *(x) = dV/dt\ x < for all x e R" if *(x) exists. 

(b) V^(x*) = if and only if x* G A, where A is the 
attractor of the dynamical system: x = f (x) . 

Here, an attractor A of a dynamical system with flow 
(fit can be fixed point, limit cycle, or chaotic attractor. In 
this sense, a potential function is a Lyapunov function. 
An attractor A is formally defined as the following. 

Definition 2 (Attractor). An attractor A of a dynamical 
system x = f (x) with flow (fit is a compact invariant set, 
with an open set U containing A such that for each x 6 U, 
(fit (x) £U for allt^O and A = f) t>0 <fi t {U). 

This definition balances between different fashions of 
literature a 20 ' 21 and is exactly the same as the definition 
of an "attracting set" in a classical textbook of dynamical 
system; 



.22 



A. Decomposition Scheme 

Many efforts generalizing Hamiltonian or gradient dy- 
namics realize that there generally exists a decomposition 
of dynamical systems x = f (x) 19 i 23 i 24 : 

x = f(x) = M(x)V*(x). 

And M(x) can be decomposed into: 

M(x) = J(x) - D(x) + Q(x), 

where J(x) and Z?(x) are semi-positive definite symmet- 
ric matrices and <3(x) is skew-symmetric. 

In our definition of the potential function ^(x), how- 
ever, matrix J(x) do not appear. That is: 

x = f (x) (2) 
= -Z>(x)V*(x) + Q(x)V*(x). 

This means that a generic system is only composed of 
an energy dissipating (gradient) part and an energy con- 
served (rotation) part. For the gradient part, potential 
^ is a common energy function; for the rotation part, ^ 
is a first integral. And we can express these two parts 
once we find the potential function for the system^: 



D = -- 



f • W 



and 



f x V* 
W • V*' 



(3) 



(4) 
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Here, / denotes identity matrix and the generalized cross 
product of two vectors defines a matrix: x x y = A = 

(^ij )nxn — {.^iUj •EjUi^nxn- 

Clearly, (BVf(x)) x V* = and (QV*(x)) • 
V\& = 0, corresponding exactly to the curl-free com- 
ponent and divergence-free component in Hclmholtz 
decomposition^. 

Further, as have been discussed in the context of 
nonequilibrium thermal dynamics^, such two parts cor- 
respond to two different structures in geometry: a dissi- 
pative bracket {•,•}; and a generalized Poisson bracket^! 

[•>•]*. 

A dissipative bracket {■, •} associates with matrix 
D(x): 

{f,g} = difDijdjg, 

and is generally defined as symmetric: {/, g} = {g,/}; 
and semi-positive definite: {/, /} ^ 0; satisfying Leibniz' 
rule: {fg,h} = f{g,h} + g{f, h}. 

While a generalized Poisson bracket [•,•] associates 
with matrix Q(x): 

[f,g] = difQijdjg, 

and can be generally defined as antisymmetric: [/, g] = 
- [g, /]; satisfying Leibniz' rule: [fg, h] = f[g, h]+g[f, h]. 

Hence, the original differential equations can be ex- 
pressed as: 

Xi = -{x h y} + [Xi,*] . 

That is, a generic dynamical system is a direct composi- 
tion of the two well-studied geometric structures. 

Later, this gradient-rotation decomposition would pro- 
vide additional insight to the understanding of chaotic 
attractors and strange attractors. 
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Figure 1 | The Simplified Geometric Lorenz 
Attractor. The dynamical system we study here is defined 
piecewisely in region Ra, Rb, and Rc, along with region 
R B i and Rc as the symmetric counterparts of Rb and Rc- 
The front, side, and top view of the regions of definition are 
shown in panel (a) through panel (c) respectively. 
Trajectories of this dynamical system would converge into 
an attractor Aj, (see subsection D and figure 7), which is a 
simplified version of the geometric Lorenz attractos^. 



by piecewise continuous ordinary differential equations 
(ODE), similar to the "geometric Lorenz attractor". We 
integrate trajectories in each continuous region of the 
model system. Then we reveal the structure of the at- 
tractor by finding the Poincarc map between the continu- 
ous regions. Through the Poincare map, the attractor is 
proved to be a chaotic attractor according to the widely 
applied definition^ of Devaney chaos. 



III. SIMPLIFIED GEOMETRIC LORENZ ATTRACTOR 

As can be seen in previous works, many efforts have 
been made to analyze Lorenz system^ as a typical model 
for chaos. Yet, to the best knowledge of the authors, 
there is only numerical evidence^ that the Lorenz equa- 
tions support a robust strange attractor—. Total un- 
derstanding of the Lorenz attractor, including but not 
limited to an analytic proof that the Lorenz attractor is 
chaotic is still lacking^!. 

An early work^ attempted to study chaotic systems by 
constructing a geometric model in a piecewise fashion to 
resemble the Lorenz system. The resultant "geometric 
Lorenz attractor" from the piecewise model is studied 
in some depth and an analogy is made between it and 
the Lorenz system^. This methodology is practically 
effective, yet the model system can become even simpler 
to be analytically proved as chaotic. 

Hence, we start out constructing a simplified geomet- 
ric Lorenz attractor. The model system is described 



A. Model System Description 

The piecewise continuous ODE model is described in 
each continuous region (from Ra to Rc, along with Rb> 
and Rc as the symmetric counterparts of Rb and Rc) 
as follows, corresponding to Figure (f ). 

f . In region R A , where x G [0,2)22, y g [-2, 2], z G [0, 2], 
yze [-2,2]: 



x = 

y = y 



(5) 



Dynamics in this region is characterized by saddle 
points at y = z = 0. These saddle points are respon- 
sible for causing bifurcation in originally close trajec- 
tories. 
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2. Region R B is defined as: x G [0 , 2/3 + 8/(3tt) x 0], 5. In region Rc, where a: G [4/3, 2], y G [-4, 1], z > 2, 
y G (2,4], z G [0,2], v^i - 2)2 + (y - 2)2 g [1,2], y/(z - 2)2 + (y + 2)2 ^ 1, ^(z - 2)2 + (y + 3/2)2 <; 

where: 5/2: 



arccos 



y-2 



V(z - 2)2 + (y - 2)2 



denoting the angle that point (y, z) form with respect 
to the center (2,2). 

In R B : 



y = 2-z 
z = y-2. 



tt/4 



(6) 



Trajectories in this region rotate for an angle of tt/2 
with respect to y = z = 2 and contract in the x direc- 
tion. 

3. In region R c , where x G [0, 2/3], y G [-1,4], z > 2, 
- 2)2 + (y - 2)2 > 1, ^(z - 2)2 + ( y _ 3/2)2 <; 

5/2: 



i = 

y = 2-z 

. 9y 21 V(3y-7) 2 + 8(z-2) 2 



(7) 



In this region, trajectories rotate for another angle of 
7r with respect to y = z = 2 and expand in the y 
direction. 

The whole system is set symmetrical with respect to 
the line: x = 1; y = 0; z G R. So we change the 
coordinate of (x, y, z) into (2 — x, — y, z) to have ex- 
pressions of the vector field in region R B i and Rc 
from expressions in region Rb and Rc- 

4. Region R B > is defined as: x G [4/3 - 8/(3tt) x 9, 2], 
y G [-4, -2), z G [0, 2], = 2) 2 + (y + 2)2 g [1, 2], 
where: 



= arccos ■ 



-y -2 



V(2-2)2 + (y + 2)2' 

denoting the angle that point (y, z) form with respect 
to the center (—2, 2). 

In R B r. 



2-x 




(8) 



9y 21 y/(3y + 7) 2 + 8(z -2p 
8 8 8 



(9) 



Vector field in region i?c"/ corresponds exactly to that 
in region Rc- 

We note that the model system in region R B ' and Rc 
is just a change of variables of the system in region R B 
and Rc- So, to avoid redundancy, we will only take 
region Ra, R B and Rc to represent all the regions of 
definition in the following analysis. 



B. Near Saddle-Focus Fixed Points 

We have constructed the model system containing one 
saddle fixed point. As in the Lorenz system, there would 
actually be another two saddle-focus fixed points when 
the system expands to the whole R 3 space. Here, we 
complete the dynamical system near the two saddle-focus 
fixed points so that the convergence behavior away from 
the attractor can be further demonstrated. 

We denote the regions of definition discussed here as 
region Rb and Rd> (see Figure [2]), each consisting of 
three parts: region -Rqa, RpB, and R D c (for region Rb 
as example). Region Rb and Rb> arc symmetrical with 
respect to the line: x = 1; y = 0; z g M, just as in 
the previous section. Hence, we follow the convention 
stated in the previous section: to take region i?/>4 , Rb& , 
and Rbc representing their symmetrical counterparts. 
The regions: Rb^, Rb b i an d Rd c an d the differential 
equations in them are written as the following. 

1. Region Rba is close to region Ra and is defined as: 
x G [0,2], y G [1,2], z G [1,2], yz G (2,4]. We simply 
take differential dynamical system in it the same as 
that in region Ra- 



x = 

z 



(10) 



-z. 



Hence, states in this region are unstable in the y direc- 
tion and stable in the z direction, causing a rotation 
effect. 

2. Region i?£>s is close to region R B and is defined as: 
x G [0, 2/3 + 8/(3 tt) x 6], y G (2,3], z G [1,2], 
^-2) 2 + (y-2)2<l. 

Here. 



Vector field in region R B i corresponds exactly to that 
in R B . 
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Figure 2 | Near Saddle-Focus Fixed Points. When the 
system expand to contain region Rd and R D > , two 
saddle-focus fixed points would emerge. This figure 
elaborates on the system near the two saddle-focus fixed 
points, (a), The regions containing the two saddle- focus 
fixed points, i.e. Rd and R D i are shown along with other 
regions, (b-d), The front, top, and side view of region Rd 
are shown respectively. 



denoting the angle that point (y, z) form with respect 
to the center (2,2). 

In region R^b : 

x 



x = - 

y = 2 - z 
z = y-2. 



tt/4 



(11) 



Same as in region Rb, trajectories in this region rotate 
for an angle of tt/2 with respect to y = z = 2 and 
contract in the x direction. 

3. In region R D c , where x <E [0, 2/3], y e [1,3], z > 2, 
A /(z-2)2 + (y-2)2 < 1: 



z + (y-2) 



z = y-2 + (z-2) 



y(z-2)2 + (y-2)2 
1 

y/( Z - 2)2 + (y 



- 1 



(12) 

In this region, trajectories tend to converge to the unit- 
radius circle centered at y = z = 2. Hence, states 
in the w hole region Rp ar e attracted to the circle: 
x = 0, ^{z-2f + {y-2f = 1. 

Here, it is observable that region Rd, containing a 
saddle-focus fixed point, form a semi-stable limit cy- 
cle at x = 0, y/(z - 2)2 + (y - 2) 2 = 1 (when y, z G 



Figure 3 Simulated Trajectory of the System. 

We simulated a trajectory of the dynamical system 
constructed. It appears to have "erratic" behaviors, similar 
to the Lorenz system. 



[1,2], the curve of the limit cycle changes expression to: 
x = 0, yz = 2). This limit cycle locates at the bound- 
aries between region Rd and its adjacent regions. This 
phenomenon corresponds with many observations that 
fixed points transit into chaotic behaviors through limit 
cycles^. 

Also, this section demonstrates that the domain of def- 
inition in the model system is not restricted to the re- 
gions discussed above. If we take dynamical system in 
the rest of M 3 space converging into the defined regions 
(Ra through Rd, Rb> through Rd'), domain of defini- 
tion can be expanded to the whole space. 



C. Trajectory and Poincare Map 

Based on Equation (5-9), we simulate the trajectory of 
the dynamical system (shown in Figure ((5)). Actually, 
trajectories in each region can be analytically solved. To 
study the structure of the attractor, we solve the trajec- 
tories in region Ra, Rb, and Rc respectively: 

1. In region Ra, trajectories are represented as: 




(13) 



where zq can usually be taken as 2. 

Hence, states in this region are exponentially unstable 
in the y direction and exponentially stable in the z 
direction. 
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The above discrete dynamical system is just a "baker's 
map" denned in literatures before^, and can be under- 
stood figuratively as the following. In the x direction, 
the mapping is contractive. The square of definition: 
[0, 2] x [—1,1] is contracted to the one third of it, form- 
ing a rectangle: [0, 2/3] x [—1, 1]. In the y direction, the 
mapping is expansive just as the doubling map22. The 
rectangle is stretched to: [0, 2/3] x [—3, 1]. Then we keep 
the right half of the resulted rectangle and move the left 
half: [0, 2/3] x [-3, -1) to the position: [0, 2/3] x [-1, 1). 
It can readily be seen that the invariant set is formed 
by iteratively removing the middle third of the intervals 
along the x direction. 



Figure 4 Poincare Map of the Dynamical System. 

This Poincare map is taken over the surface of z = 2. It can 
readily be observed that this map is a "baker's map"— and 
creates a Cantor set multiplying a real line segment. 



2. In region i?p, trajectories are: 



1 4 



V = \ /vl + z l cost + 2 
z = "Jvl + Z l shit + 2, 



(14) 



where yo can be 2. 

Here, we can observe that x(t) decreases monotoni- 
cally while y(t) and z(t) form a circle. 

3. In region Rc, trajectories are: 

X = Xq 

V = Vvi + 4 cost + - (l - y/yj + z 2 ) + 2 (15) 

z = VVO + z o sint + 2 : 

where zq can be 2. 

Hence, trajectories in region Rc move along circles 
determined by initial conditions. 

Then, to further study the structure of the attractor 
of the system, we want to calculate the Poincare map of 
the system. Here, we take Poincare surface of section as: 
z = 2, and find the resultant Poincare map as a discrete 
dynamical system defined on [0,2] x [—1,1] (shown in 
Figure (jlj and follows). 

When (x,y) G [0,2] x [0,1], 



3 Xr 



r i 

J x n+1 = 

{ Un+1 = 2?/„ 

When (x,y) 6 [0,2] x [-1,0), 

1 




(16) 



(17) 



D. Attractor of the Model System 

We denote the attractor of the model system as: Aj,. 
And with the Poincare map of the model system defined 
as a dynamical system on [0, 2] x [—1,1] (in the previous 
section), we denote its attractor as: Ap. In this sec- 
tion, we first express attractor Ap of the Poincare map 
in terms of the Cantor set; then we express attractor A^ 
of the original model system. 

We have already found the Poincare map of the model 
system as iteratively removing the middle third of the in- 
variant sets along the x direction. That is: first, remove 
the set (2/3, 4/3) x [—1,1]; then, remove the middle third 
of the left two sets [0,2/3] x [-1, 1] and [4/3,2] x [-1, 1]; 
and iterate the process all along. We present all the re- 
moved intervals iteratively as the following: 



I I.TT ) x [-1,1] , 



3' 3 



and 



y^ + 4^ 



Then the attractor of the Poincare map is [0, 2] x [— 1, 1] 
minus the union of all the sets Cf. 

oo 

Ap = [0,2] x [-1,1]- IjQ x [-1,1] 



[0,2]-(jQ x [-1,1] 
v i=l / 

:x [-1,1], 



(19) 



where C denotes the Cantor set^ defined on the interval 
of [0,2]. 

Hence, attractor Ap of the Poincare map is the Cantor 
set multiplying a real line segment. We can calculate its 
box-counting dimension^ to be: 

4(A P ) = lim 1o & N ^ A p) = ! + j n ( 2 )/Zn(3). (20) 
e^o log(l/e) 



7 



Hence, the attractor Ap is of fractal dimension, a strange 
attractoi— . 

With the trajectories of the system analytically 
solved in each region, we further express attrac- 
tor Al of the model system as (assuming 9 = 
arccos (y - 2)/y/(z - 2) 2 + (y - 2) 2 ): 

In region Ra, ifC; 

in region R B , (4/(3tt) x 6 + 1/3) -1 x E C; 
in region Rc, ieC. 

The box-counting dimension of attractor Al is then 
calculated to be: 



lim 



logjV(e,A £ ) 
log(l/e) 



= 2 + ln(2)/ln(3). (21) 



It is hence a strange attractor with fractal dimension. 

It will be proved in the following section that attractors 
Ap and Al are also chaotic attractors. 



E. Proof of the Attractor as Chaotic 

By the widely applied definition^ of Devaney chaos, 
an attractor A is defined as a chaotic attractor if: 

1. the attractor is indecomposable (i.e., if ^ A' C A is 
an attractor, then A' = A); 

2. the system is sensitive to initial conditions when re- 
stricted to A (defined in the following definition [3]) . 

Definition 3. A map (a continuous-time system is de- 
fined similarly) has sensitive dependence on initial con- 
ditions when restricted to its invariant set A, if there 
exists r, for any po G A, and S > 0, there exists p' Q G A: 
\p'a — Po I < an d an iterate k > such that 



\f k (Po)-f k (po)\^r. 



(22) 



Attractor Al has already been taken as the smallest 
attracting set, so it is an indecomposable attractor by 
default. We just need to prove that the system has sensi- 
tive dependence on initial conditions when restricted to 
A L . 

We first prove that attractor Ap of the Poincare map 
is chaotic using the fact that the doubling map22 is sensi- 
tively dependent upon initial conditions when restricted 
to its attractor. Then we prove in exactly the same way 
that attractor Al of the model system is chaotic by the 
sensitivity of Ap. 

Proposition 1 (Sensitive Dependence of the Poincare 
Map). The Poincare map of the model system has sen- 
sitive dependence upon initial conditions when restricted 
to its attractor Ap. 

Proof. For any po G Ap , with its neighboring initial point 
p' G Ap, we take p' as (x' ,y' ) = (x' ,y' - 5 ■ Sign(y )). 



Clearly, \\p n - p'J > 
to initial conditions of 
Doubling Map: 



Vn -y' n \- So, proving sensitivity 
Ap is equivalent to that of the 



Z/n+l = L2j/nJ, (Vi G (0,1), Vi). 



(23) 



With the sensitive dependence upon initial conditions of 
doubling map when restricted to its attractor proved^, 
Poincare map of the model system is also sensitive when 
restricted to its attractor Ap. □ 

In exactly the same way, the model system can thus 
be proved sensitively dependent upon initial conditions 
when restricted to its attractor A^. Hence, it is a chaotic 
attractor by definition. 

We also calculate the commonly used indicator of 
chaos: Lyapunov exponents^ for the model system at 
fixed points. By solving the Lyapunov exponents in each 
coordinate direction, we find that in region Ra- l x = 0, 
l y = 1, and l z = — 1. It can be found that there is a 
positive Lyapunov exponent l y = 1 denoting exponen- 
tial expansion in the y direction. In region Rb and Rq, 
i x = l y = l z = 0, which means that the expansion effect 
causing the sensitivity of the system is mainly exerted in 
region Ra- 

From this and the former section, we find that the 
attractor of the model system is a chaotic attractor with 
fractal dimension: a strange chaotic attractor—. 



IV. CONSTRUCTION OF POTENTIAL FUNCTION IN 
THE CHAOTIC SYSTEM 

Based on the above observation, we want to start 
constructing a potential function to describe the over- 
all structure of the chaotic system. First, we construct 
a "seed function", J 7 , to account for the "strangeness" 
of the system's attractor. Then we prove its continuous 
differentiability so that it can be applied in the construc- 
tion of potential function for the model system. Later, 
we explicitly express the potential function in terms of 
the seed function T . 



A. Definition of the "Seed Function" T 



Definition 4 (Function T). Let 



0, xG [0, 2/3] (J [4/3, 2] 
1-cos(3ttx), x€ (2/3, 4/3); 



and 



l/9x/ n (3x), xe [0,2/3] 
fn+i{x) = { 0, xe (2/3, 4/3) 

1/9 x f n (3x-4) , x G [4/3, 2] 

Thus, we define the function T(x) as: 

oo 



(24) 



(25) 



(26) 
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Figure 5 Illustration of the Self-Similar Function: 

T. We hereby plot function T to intuitively visualize it. 
With T, construction of a potential function in the chaotic 
system would be natural. 



Function J-(x) defined on [0, 2] is shown in Figure (|5|). 
It has a fractal structure as the attractor Ap of the 
Poincare map. 



B. Proof of T(x) as Continuous Differentiable 

Here, we propose that the function J-(x) defined above 
is continuously differentiable. 

Proposition 2 (Continuous Differentiability of Function 
J-). Function J-(x) = X)n=i/"( x ) defined in Definition 
2] is continuously differentiable. 

Proof. We first calculate the derivative of every /„ (x) on 
[0,2]. Then we bound f n {x) and \f' n {x)\ by geometric 
series to prove uniform convergence of their sums, and 
hence prove that F(x) = Y^Li fn{x) is continuously dif- 
ferentiable. 

From the definition (in equation (|24p and ([2"5]0 : 




0, xe [0, 2/3] U [4/3, 2] 

1 -cos(3ttx), x G (2/3, 4/3) . 



l/9x/„(3x), xe [0,2/3] 

0, x £(2/3,4/3) 

1/9 x /„(3x - 4) , a; e [4/3, 2] 



we have by taking derivative on both sides: 
fi(x) - 



0, xG [0, 2/3] (J [4/3, 2] 
3tt sin(37rx) , x € (2/3 , 4/3) , 



and 



l/3x#(3x), xe [0,2/3] 
0, xG (2/3, 4/3) 
l/3x/;(3i-4), a; G [4/3, 2] 



Obviously, fi(x) < 2 and [/{(x)) < 37r. And also, 
fn(x) is continuous for any x G [0, 2]. 

If we further denote the set in which f n {x) is nonzero 



we have: 



2 4 
3' 3 



and 



Cn+l 



(27) 



Since C„ C [0,2], 



c 



3 1 1 3 

We thus can conclude that 
1 



n 



«S g/ B (x) ^ 9-"/i(^) < 2 x <T", 



and 



< < 3-"|/{(a;)| < 3^ x 3"", 

for any x G [0, 2]. 

At this point, we give an upper bound for the series 
Y^=i fni x ) an d its derivative (Although the least up- 
per bound is even smaller if we note that <C n f| C n +i is 
actually empty, an upper bound is good enough): 



and that 



m Q 1 



Ei£(*)i<£-T* 3 - 



So, ~Y^=if' n {x) is uniformly absolutely-convergent. 
Hence, J2n=ifn( x ) and X)£Li fn{x) are all uniformly 
convergent. 

With every f n {x) continuously differentiable, J~(x) = 
Sn=i fn(x) is continuously differentiable: 



rt_ 

d.r 



(28) 



□ 



Clearly, the points where F(x) = J~'(x) = form a 
Cantor set C corresponding to the attractor A^ of the 
model system (J-(x) = if and only if a; G C). At this 
point, we found that the potential function can be con- 
structed in the following section. 
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C. Constructing Potential Function in the Chaotic System 

As in equation (fTTj) we use 9 to denote the angle that 
(y,z) form with respect to the center (2,2): 



9 = arccos 



tf-2 



Then, we construct potential function in each region re- 
spectively: 

1. In the right part of region Ra, where x,y,z G [0,2], 
yze [-2,2]: 



(29) 




2. In region R B , where x G [0,2], y G (2,4], z G [0,2], 
V(z-2) 2 + (j/-2)2e[l,2]: 

3. In region where x G [0,2/3], y G [-2,4], z > 2, 
V(z-2) 2 + (y-2)2 > 1, v/(z - 2)2 + (y - 3/2)2 <; 
5/2: 

• -(-s» + S)^»- (31) 

Here, we plot the potential function taken on the 
Poincare section in Figure (|6]). 

If we are also interested in the dynamics near saddle- 
focus fixed points, potential function in region R^a, 
Rd b i and Rjjc can also be constructed as follows. 

1. In region R d a, where x G [0,2], y G [1,2], z G [1,2], 
yze [2,4]: 

^ = {l e+ l)^ + 1 -{l yZ - 2 ) 2 - (32) 

2. In region R d b , where a; G [0, 2/3 + 8/(3tt) x 9], y G 
[2,3], zG [1,2], ^-2) 2 + (y-2)2 e [0,1]: 

+ 1 - (z I f - (y 2) 2 . 

3. In region R^a, where x G [0, 2/3], y G [1,3], z > 2, 
V^i - 2)2 + (y - 2)2 <; 1: 

^=(-±9+^(3*) (34) 
+ 1 - [z - 2) 2 - (y 2) 2 . 



Figure 6 Potential Function on Poincare Section. 

On the Poincare section, the potential function is 
demonstrated to be a fractal object: it is zero when point 
(x,y,2) belongs to the attractor Al, i.e., \1/|^ = 2 = when 
ieC. And when (x,y, 2) does not belong to the attractor, 
the potential function ^ has self-similar structure. 



Potential function in region Rp is gradually higher in the 
center of the region than in its boundary with other re- 
gions. Hence, points in this region will naturally converge 
to its boundary (at: x = 0, y/(z - 2) 2 + (y - 2) 2 = 1, 
when y > 2 or z > 2; and at: x = 0, yz = 2, when 

y,ze[i,2}). 

The construction of the potential function is not 
unique. If we change the expression of fi(x) in the defini- 
tion of the "seed function" J-(x), we can have a different 
potential function for the system. 



V. VERIFICATION OF THE POTENTIAL FUNCTION 

Here, in this section, we verify the integrity of the po- 
tential function from three angles: First, we show that 
it is continuous in the domain. Second, we demonstrate 
that it decreases monotonically along the vector field. 
Third, we show that V^(x) = if and only if x belongs 
to the attractor of the system. 



A. Continuity of the Potential Function 

1. At the boundary of region Ra and Rb, y = 2, x G 
[0,2], zG [0,1]. 

Hence, in equation (|29l) . 9 = arccosO = 7r/2; 
(4/(3tt) x 6 + 1/3)" 1 = 1; and (4/(9tt) x 9 + 5/9) = 



10 



7/9. So, 



Since x = and y, z G [0, 2], 



*A|«=a- = g-^O) = *s 1 | !/ =2+- 

2. At the boundary of region i?s and Rc, z = 2, x G 
[0,2/3], and j/ £ [3,4]. 

Hence, in equation (f3"U)l , = arccosl = 0; 
(4/(3?r) x 6» + 1/3)" 1 = 3; and (±4/(9tt) x + 5/9) = 
5/9. So, 

3. At the boundary of region Rc and Ra, z = 2, x G 
[0,2/3], and y G [-1,1]. 

Hence, in equation (|3"Tj) , = arccos( — 1) = 7r: 
(4/(9vr) x 9 + 5/9) = 1; and (-4/(9vr) x 9 + 5/9) = 
1/9. So, 

*aU=2" 

It can be observed that in the definition of J-(x) = 

fn+l(x) = lfn(3x), XE [0,2/3]. 



Clearly, because /i(a;) = in [0,2/3], 

-7-(3x) = ^/ n (x)=J-(x), 



(35) 



n — 2 



for a; G [0,2/3] 



Please note that this is a critical point in the construc- 
tion of potential function in cases of self-similar attrac- 
tors. If the potential function constructed is not self- 
similar accordingly, then the boundary would not totally 
fit. 



4 (2 - y)z + (2 - z)y 



F{x) < 0. (36) 



9tt (z - 2) 2 + (y- 2) 2 ' 

2. In region R B , where x G [0,2], y G (2,4], z G [0,2], 
v/( Z -2)H(ff-2)'6[l,2] : 



/ y-2 
:r = — x — h arccos — 

1 4 v/FWIFf 



7T 



z-2 



r/7 



y-2 



(z-2)2 + (y-2)2 a (z-2)2 + (y-2)2 
-1. 



9tt 



3tt 3 



1 



4 1\ 4 • 
3^+3 ^3^ 



4 

3^ 



9t+ 



-J" 



4 1 

3^ + 3 



±0 + ^ fA, 

9tt 9/ V3tt 
3?r 3/ U 



4 

37' 



-i 4 
x H x 



3tt 



Since 



±0+i) r- 

3tt 3 / V4 



1 4 

IE ■ 



3tt 



x = 0, 



9?r 



*b = ~z — J~ 11—0 + -) x \ < 



3tt 



(37) 



B. Monotonic Decreasing of the Potential Function 

Then we take Lie derivative (derivative along the vec- 
tor field)22 of the potential function along vector fields in 
each region remembering that J-(x) ^ for any x G [0, 2]. 

1. In the right part of region Ra, where x,y, z G [0,2], 
yze [-2,2]: 



( x = 

6 = Sgn(z - 2) 



< - 



z-2 



y-2 



(2-2)2 + ^-2)2^ (z -2)2 + (y- 2)2' 

(2 - y)z + (2 - z)y 
(z - 2)2 + (y - 2)2 • 



3. In region R c , where x G [0,2/3], y G [—1,4], z > 2, 
V(* - 2) 2 + (y - 2)2 ^ 1, ^(z - 2)2 + (y - 3/2)2 ^ 
5/2: 



r ± = o 

= - 



z-2 



y 2 



(2-2)2 + ^-2)2^ (z -2)2 + (y- 2)2- 

(z-2) 2 y-2 



( (z - 2)2 + (y -2)2 (z-2)2 + (y - 2)2 
9y 21 ^(3y-7) 2 +8(z-2) 2 



> 0. 



Hence, 



J"(3x)0 < 0. 

9tt 



(38) 
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C. Potential Function and the Attractor 



1/2! € [-2,2]: 



Here, wc verify that the potential function attains ex- 
tremum: V v E'(x) = if and only if x = (x, y, z) be- 
longs to the att ractor A^ of the s ystem. Again, 9 = 
arccos (y - 2)/y/{z - 2) 2 + (y - 2) 2 . 



4 z-2 

9tt (z - 2) 2 + (y-2) 
4 y-2 

\ 9tt (2 - 2) 2 + (y - 2)2 U 

Here, .F(x) = and J-'{x) = if and only if x 6 C; and 
point (x, y, 2) belongs to the attractor A^ if and only 



jjuim yju^y^) LfUlUll^S LU L11U clLLl cHJLUI 11 CL11L1 u 

1. In the right part of region Ra, where x,y, z G [0, 2], if x € C. So, V^a = if and only if (x, y, z) G A^ 



2. In region R B , where x G [0,2], y G (2,4], z G [0,2], ^(z - 2) 2 + (y - 2) 2 G [1,2]: 




^ 9tt (z - 2) 2 + (y - 2) 



/ 



3. In region C, where x G [0,2/3], 2/ G [-1,4], z > 2, 
v /(z-2) 2 + (y-2) 2 > 1, ^(i - 2) 2 + (y - 3/2) 2 < 
5/2: 



Wf7 = 



■A* 

3tt 



z-2 



\ 



9tt (j? - 2)2 + (y - 2)2 
4 y-2 



\ 9tt (z - 2) 2 + (y - 2) 



J"(3x) 
,-.F(3x) 



Here, J" ((4/(3?r) x 9 + 1/3) 1 xj = and T' ((4/(3tt) x 9 + 1/3) 1 xj = if and only if (4/(3tt) x 9 + 1/3) X xG 

C; and point (x, y, z) belongs to the attractor Al if and only if (4/(37r) x 9 + 1/3) 1 x G C. So, Vf b = if and 
only if (x, y, z) G A/,. 

I 

approach^, energy-like function technique^, minimum 
action method^, and etc. These previous methods all 
construct a potential-like scaler function to analyze cer- 
tain chaotic system. Unfortunately, the scalar functions 
in these works all lack certain important properties. 

For example, the generalized Hamiltonian systems ap- 
proach takes a quadratic form of the state variables as 
the "generalized Hamiltonian"—, a Hamiltonian that in- 
cludes conserved dynamics, energy dissipation, and en- 
ergy input. The third part transforms an autonomous 
differential equation into a non-autonomous physical 
model, attempting to explain for the "irregular"— be- 
havior of chaotic systems. Contrary to this expectation, 
when the Hamiltonian is set possible to dissipate and 
increase, the generalized Hamiltonian itself becomes a 
chaotic oscillating signal with respect to time. There- 
fore, it remains an issue as to what additional insight this 
generalized Hamiltonian can provide about the original 
system, such as global stability or local performance. 

The energy-like function technique is essentially simi- 
lar to the generalized Hamiltonian approach. Its energy- 
like function differs from the generalized Hamiltonian in 
a way that it may not be a quadratic form of the state 
variables. Rather, the energy-like function is constructed 
based on the "geometric appearance"— of the attractor 
corresponding to the specific chaotic system. Although 
this technique would seem more sophisticated, its encrgy- 



Hcrc, .F(3x) = and .F'(3x) = if and only if x G C: 
and point (x,y,z) belongs to the attractor A^ if and 
only if x G C. So, V^c = if and only if (x, y, z) G 
A L . 

In exactly the same way, can we also show the integrity 
of the potential function in region Hp. 



VI. COMPARISON WITH RELATED WORKS 

As discussed in the introduction, constructing a 
potential-like function in the chaotic system is actually an 
effort that is by no means totally strange to researchers. 
Until recently, there are various efforts seeking to de- 
scribe chaotic dynamics using generalized Hamiltonian 



12 




(a) Front View 


(b) Side 


(c) Top View 


View 



Figure 7 | Strange Chaotic Attractor. 

We find that a connected surface of the attractor is of 
infinite layers. We show how surface x = 2 is linked to the 
other equipotential layers. It can be seen that the surface of 
the attractor is orientable. We hereby demonstrate the 
strange chaotic attractor viewed from (a), front; (b), side; 
and (c), top. The trajectory running from point (2, 1/4, 2) 
is also shown in the figure. 



like function still oscillate chaotically along with time, 
describing chaotic dynamics in a chaotic fashion. Loss 
of monotonicity restricts the function from describing 
the system's essential properties like stability and per- 
formance. 



The minimum action method, however, cast the prob- 
lem under the light of zero noise limit. By construct- 
ing an auxiliary Hamiltonian^ (commonly denoted as 
"Frcidlin-Wcntzcll Hamiltonian"), Frcidlin-Wentzell ac- 
tion functional can be minimized^. This method ana- 
lyzes chaotic system by possible transitions between limit 
setsi£. But since the Frcidlin-Wcntzcll Hamiltonian can 
be not bounded even in globally stable systems, it is not a 
quantitative measure comparable between points in state 
space, hence, not an ideal potential function. 

In short, all the previous works each focuses on one 
attribute of the potential function. However, as we can 
see from our constructive result, only when all the re- 
quirements (in definition 1) are met, would the potential 
function reflect evolution of the whole system and struc- 
ture of the chaotic attractor. In this sense, the current 
work is the first construction to satisfy such strong con- 
ditions, providing a both detailed and global description 
for a chaotic system. 



VII. CHAOTIC ATTRACTOR AND STRANGE 
ATTRACTOR 

With the potential function constructed, we can easily 
solve the system's attractor without any need of numeri- 
cal simulation. Wc find that the attractor is composed of 
connected surfaces, each of infinite layers. Starting from 
the plane x — 2, we show the configuration of these layers 
in Figure ([7]) . Since geometric configuration of chaotic at- 
tractor interests many researchers^, we demonstrate in 
the figure that the chaotic attractor of the system studied 
in this paper consists of orientable surfaces. 

The chaotic attractor here is a strange attractor of frac- 
tal dimension^. And in the literature of dynamical sys- 
tems, there have long been discussions about the relation- 
ship between chaotic attractors and strange attractors^k 
Several examples of strange nonchaotic attractors and 
nonstrange chaotic attractors have been found^. Until 
recently, strange nonchaotic attractors are still studied^. 

The potential function approach provides a unified 
framework to treat the topics of chaotic attractors and 
strange attractors together. To clarify this insight, we 
need to apply our decomposition method (equation ([5])) 
here: 



x = f(x) = -DVtf (x) + QVf(x). 



where 



and 



D 



f ■ W 



f x V* 
^ ~ W • V*' 

We first analyze our model system with this decompo- 
sition framework. Then we further modify our model sys- 
tem to two typical cases interesting to many researchers: 
a nonstrange chaotic attractor and a strange nonchaotic 
attractor. After analyzing these two cases, we explain 
the different origins of chaotic attractors and strange at- 
tractors in general. 



A. Decomposition of the Chaotic System 

According to our decomposition scheme, we first de- 
compose the chaotic dynamical system in each region 
into two components: the gradient component and the 
rotation component. 

In region Ra, V'J'a is solved as: 



I (—1 

\9ir 
4 



\ 



9?r (z 



2) 2 
2 - 



V 9tt (i 



(y 
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Hence, we can find the expression of the matrix Da 
accounting for the gradient component of the vector field 

in region R A : D A = — —I (39) 



4 {2 - y)z + (2 - z)y 
9tt (z-2) 2 + (y-2) 2 



The decomposed gradient part would then be: 

4 (2 - y)z + (2 - z)y 



D A V^ A 



9tt - 2) 2 + (y- 2) : 



.F(x) 



4 2-z 



9tt 



(J-'(x)) 2 



(4/9tt)' 



(z - 2)2 + (y - 2)- 



9tt (z-2) 2 + (y-2) 2 ^ 

V 9^(z-2)2 + (y-2)2^ W / 



Also, we can find the decomposed rotation part by finding Qa &s: 
f a x 



1 



9tt 9 
4 



(z-2 ( )2 + (l-2)2 (J " (a;)r 




V 



o \ . . . 4 (2 - z)z — (2 — y)y . 

,9tt 9/ W 9tt (z - 2)2 + (y - 2)2 v ; 



4 (2 - z)z - (2 - y)y 

97T (2 - 2)2 + (j, - 2)2 





.F(x) 



The decomposed rotation part would then be: 



Qa^^a = 



1 



4 ^ 5 Vr7r// (4/9^) 2 

^ + 9j + (z-2)2 + (,-2)2 



(.F(z)y 



/ 4 



^ fa-2> + (^-2) V ( ^ (a;) 



9tt V9tt 9 J (z - 2) 2 + (y - 2) 2 
4. 5\ 2 ,_, 2 (y-2) 2 y-(y-2)( 2 -2> 
9tt 



i^n) (^)) 2 



((*- 2)2 + (1,-2)2)' 



(^(a:))' 



9^+9 <^<*» + U?J 



4 V (y-2)y(z-2)-(z-2) 2 z 



((z-2)2 + (y-2)2r 



(40) 



(41) 



When the system approaches its attractor, i.e., F(x) —> 0, 



-] (l-cos(37rx)) 

JCW = lim _W (42) 

(^'(x)) 2 ffl y s 

— J 07rsm(37rx) 

1\" 



lim 



(3ttx) 2 
97 2 V ; 



x-To f f\\ n \ 2 18tt 2 ' 

97T 2 X 
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Thus, when the system converges to its attractor, the 
gradient matrix Da would be: 



D A 



4 {y-2)z+{z~2)y 
9tt (z - 2) 2 + (y - 2) 2 .F(x) 
_4_ 5\ 2 (-F'(x)r 

(y - 2)z + (z - 2)y 

(z - 2) 2 + (y - 2) 2 r 



(43) 



81tH 



9tt 



which is finite. 

Hence, the gradient component DaV^a of the system 
would converges to zero when approaching the attractor. 
So the motion on the attractor is caused totally by the 
rotation part: Qa^^a- 

In exactly the same way, the decomposition procedure 
can be carried out in region B and region C, and the same 
conclusion holds. 



B. Nonstrange Chaotic Attractor 

Let's first examine an example of nonstrange chaotic 
attractor by modifying our original system a little (in 
region R B , equation pTT) ): 

In region Rb, we set 

9 = arccos 



^/(z-2) 2 + (y-2) 2 ' 



as in equation (fTTj) . Then we change the dynamical sys- 
tem in region Rb (defined as x G [0 , 4/it x 6], y € [2,4], 
z G [0, 2], y/{z - 2) 2 + (y - 2) 2 G [1, 2]) to: 



x 

x = — — 



y 

z ■ 



The same as in the original system, domain of definition 
can be expanded to the whole K 3 space. 

Consequently, the Poincarc map would be as follows: 

When (x,y) G [0,2] x [0,1], 



X n +1 = 

y n+ i = 2y n - 1. 



When (x,y) G [0,2] x [-1,0), 



X n +1 = f 

y n +i = 2y„ 



The attractor of the modified system would be: 
(assuming 9 = arccos (y — 2)/y / (z — 2) 2 + (y — 2) 2 ): 
In region Ra, x — or 2; 
in region R B , (tt/2) x (x/9) = or 2; 



in region Rq, x — or 2. 

The attractor is shown in Figure ([5]) . We can calculate 
its box-counting dimension to be: 



d b (A' L ) = lim 



logJV(e,Ai 



log(l/e) 



= 2, 



(44) 



which is an integer dimension. Actually, the attractor A' L 
is just two orientable surfaces folded together. Hence, it 
is no longer a strange attractor anymore. 

Exactly as in the original system, the modified attrac- 
tor can be proved to be chaotic. And we can also cal- 
culate the commonly used indicator of chaos: Lyapunov 
exponents 2 ^ for the model system at fixed points. Lya- 
punov exponents are solved in each direction as: l x = 0, 
ty = 1, and l z = — 1 in region Ra (In other regions, 
i x = £ y = l z = 0). It is found that there is a positive 
Lyapunov exponent £ y — 1 denoting exponential expan- 
sion in the y direction, exactly as in the original model 
system. 

So, it is clear that the modified attractor is a non- 
strange chaotic attractor. 

Now, we construct a potential function $ for the new 
dynamical system by first appointing a new seed function 
F(x) defined in [0,2]: 

F(x) = 1 - cos(tt:z;) , x G [0 , 2] . 

Hence, the potential function can be represented as: 

1. In the right part of region Ra-, where x,y,z G [0,2]: 



F(x). 



2. In region R B , where y G [2,4], z G [0,2], 
V(^-2) 2 + (y-2) 2 G [1, 2], x G [0 , 4/tt x 9}: 



3. In region Rc, where x = 2, y G [—2,4], z > 2, 
y/(z - 2) 2 + (y - 2) 2 ^ 1, v /(z-2) 2 + (y-3/2) 2 < 
5/2: 

*c = 0. 



Here, $ = corresponds to the attractor. 
We can further decompose the system as with the orig- 
inal model system: 

x = f(x) = -L>V$(x) + QV$(x). 

Then, in region Ra'- 



/ 9F'(3 



V$4 = - 



2 - z 



F(x) 



(z-2) 2 + (y-2) 2 

. Fix] 

\ (z-2) 2 + (y-2) 2 ^ W 
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D 



(2-y)z+(2-z)y 
(z - 2)2 + (y - 2)2 



(45) 



2 (-F'(x))^ 



1 



(F(x)Y 



(z - 2)2 + (y - 2)2 
The decomposed gradient part would then be: 
(2-y)z+(2-z)y 



(z-2)a + (y-2): 



-F(x) 



6 2 (F'(x)Y 
/ 6F'{x) 



(F(x)Y 



(z-2)2 + (y-2)a 



2 - z 



(2 - 2)2 + (y - 2)2 

. IzJL F ( X ) 

V (z-2)2 + (y-2)2^ W 



(46) 



J 



Also, we can find the decomposed rotation part by finding Qa as: 



M = 



2 (F'(x)) 2 + (F{x)) 2 



( -yOF'(x) 
y9F'(x) 

(2 - z)z - (2 - y)y 



\ 



-z6F'{x) 



(2-2)2 + ^-2)2 



F(x) 



z6F'{x) 
(2 - z)z - (2 - y)y 
(z - 2)2 + (y- 2)2 





The decomposed rotation part would then be: 
QaV^a 



(47) 



^ (F'(x)r+(F(x)y 

( (y-2)z + {z-2)y aT?u 
' - (,-2)2 + ^-2)2 ^ 
^2(^^)2., ( y -2) 2 y -(,-2 )( ,-2)z 



(48) 



(( 2 -2)2 + (y-2)2)^ 



(F(x)Y 



-z9 2 (F'(x)) 



2 , (y-2)y(z-2)-{z-2Yz 
(( 2 -2)2 + ( y -2)2) 2 

I 



(F(x))^ 



The properties of the gradient part and the rotation We simply take the gradient of potential function of 
part corresponds exactly to the original model system. the original system in each region of definition: 
That is: the gradient part converges to zero when ap- 
proaching the attractor; motion on the attractor is de- 
termined by the rotation part. ( ^ _ _g ^ 

y = -d y V 
z = -d z ^. 

C. Strange Nonchaotic Attractor 



A strange nonchaotic attractor can also be con- 
structed. 



If we take left part of region Ra for example, the vector 
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(a) Front View 


(b) Side 


(e) Top View 


View 



Figure 8 | Nonstrange Chaotic Attractor. 

We change the expression of the system a little, so that the 
attractor is just two orientable surfaces folded together, 
rather than a fractal structure. However, it remains to be a 
chaotic attractor. We hereby demonstrate the nonstrange 
chaotic attractor viewed from (a), front; (b), side; and (c), 
top. The trajectory running from point (2, 1/4, 2) is also 
shown in the figure. 



field would be: 




9^ (z-2) 2 + (2/-2) 2J " ( - T) 
o 



\ 9tt 



(y 



J 



The resultant ODE system defined by the gradient is a 
dynamical since it is Lipschitz continuous in each region. 
And with existence and uniqueness of the flow guaran- 
teed by Lipschitz continuity, conditions for the system 
being a dynamical system can be satisfied and extended 
to include boundaries. 

The system would converge downward the potential 
function ty until reaching the states where ty = 0. Con- 
sequently, attractor of this system is characterized by 
^ = 0, as in the original model system. Hence, the gra- 
dient system's attractor is the same attractor Al of the 
original model system, whose box-counting dimension: 



logA(e,A £ ) 
d b (A L ) = hm — ■ — — — — = 2 
e^o log(l/e) 



ln(2)/ln(3). (49) 



Hence, the system has a strange attractor. 

Since = when VP = 0, the dynamical system is 
not sensitively dependent upon initial conditions when 
restricted to the attractor. So, the attractor is not 
chaotic. Also, its Lyapunov exponents at the fixed points 
(where ieC,i/ = z = 0) would be: 4 = -(V 3 )" x 97I " 2 > 
t y — i z = 0. Hence, it's a strange nonchaotic attractor. 



Decomposition of this system would give: D(pt) = I 
and Q(x) = 0. Thus, DV*(x) = -f(x) is just the re- 
versed gradient system. 



D. Chaotic Attractor versus Strange Attractor 

The previous two examples show that the concepts of 
chaotic attractor and strange attractor do not imply each 
other. Under our framework of decomposition (equation 
©) here: 

x = f (x) = -DVtf (x) + QV*(x). 

Since 4> = V* • x = V* • (DVf), * decreases mono- 
tonically according to the gradient component DVty of 
the vector field /. Then the attractor is naturally char- 
acterized by DVf = 0. So, whether the attractor is a 
strange attractor is determined by the gradient part of 
the vector field. 

Sensitive dependence upon initial conditions when re- 
stricted to the attractor, however, is determined by the 
rotation part of the vector field: QV\P. Once the sys- 
tem has evolved to the limit set, DVty would equal to 
zero, and QVty would be prevalent. Hence, the rota- 
tional vector field on the attractor causes the expansion 
of the state space, leading to dynamical sensitivity. Con- 
versely, when Q\7^> = 0, chaotic motion on the attractor 
would not exist. In this sense, nonzero rotation part of 
the dynamical system is a necessary condition for causing 
hyperbolic chaos^i. 

To sum up, gradient part and rotation part of the vec- 
tor field are responsible for the creation of strange at- 
tractor and chaotic attractor correspondingly. Although 
they are both affected by the geometrical configuration 
of the potential function vp, they denote dissipation and 
circulation respectively. 



VIII. CONCLUSION 

In the present paper, it is shown that potential func- 
tions with monotonic properties can be constructed in 
continuous dissipative chaotic systems with strange at- 
tractors. The potential function here is a continuous 
function in phase space, monotonically decreasing with 
time and remains constant if and only if limit set is 
reached. This definition is a natural restriction of generic 
dynamics since it is a direct generalization of Lyapunov 
function and corresponds to the concept of energy. 

Potential function defined this way also implies that 
the dynamics can be decomposed into two parts: a gra- 
dient part, dissipating energy potential; and a rotation 
part, conserving energy potential. The gradient part 
drives the system towards the attractor while the rota- 
tion part perpetuates the system's circular motion on the 
attractor. 

To demonstrate the power of this framework in chaotic 
systems, we simplify the geometric Lorenz attractor, and 
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prove by definition that it is a chaotic attractor. Then we 
analytically and explicitly construct a suitable potential 
function for the attractor, which, to our best knowledge, 
is the first example in chaotic dynamics. The potential 
function reveals the fractal nature of the chaotic strange 
attractor. 

Wc further analyze the concept of chaotic attractor and 
strange attractor with our decomposition. It is found 
that chaotic attractor originates in the rotation part, 
prompting the state space of the attractor to expand; 
while strange attractor originates in the gradient part, 
causing initial states attracted to complex limit set. 
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